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and  total  eddy  heat  transports,  computed  as  functions  of  latitude  and  long¬ 
itude  for  three  individual  Januaries,  are  described  and  discussed.  Corre¬ 
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CHAPTER  ONE 
INTRODUCTION 


1.1  Background 

In  the  troposphere,  the  equator  is  warmer  than  the  poles,  and  this 
is  a  result  of  radiational  imbalances .  We  observe  a  net  heating  due  to 
radiation  at  the  equator,  and  at  the  poles  a  net  radiational  cooling  occurs. 
If  this  was  the  only  process  working,  the  equator  to  pole  temperature  gra¬ 
dient  would  change  and  become  larger.  Looking  at  a  long-term  average,  the 
temperature  structure  of  the  atmosphere  is  essentially  in  equilibrium.  For 
this  condition  to  exist,  heat  must  be  transported  poleward. 

In  the  atmosphere,  the  south-to-north  transport  of  sensible  heat 
at  a  given  time  and  location  can  be  expressed  mathematically  as 

HT  =  (pC  AZ)VT  (1.1) 

P 

where  HT  is  the  local  south-to-north  sensible  heat  transport,  V  is  the 

south-to-north  wind  component,  T  is  the  absolute  temperature,  C  is  the 

P 

specific  heat  of  dry  air  at  constant  pressure,  p  is  the  density  of  the 
ear,  and  AZ  is  the  thickness  of  the  layer  in  which  the  heat  transport  is 
to  be  computed.  In  developing  the  thermodynamic  equation  for  practical 
use  and  study,  this  transport  is  often  averaged  over  an  appropriate  time 
interval  and  averaged  zonally  (i.e.  averaged  around  a  latitude  circle). 

This  averaging  may  be  represented  symbolically  by: 

[HT]  =  (pC  AZ)  [VT]  (1.2) 

P 


where  the  bar  represents  a  time  average  and  the  brackets  represent  an 
average  around  a  latitude  circle.  Equation  1.2  represents  the  mean  meri¬ 
dional  heat  transport.  This  transport  is  commonly  broken  down  into  three 
components  by  expressing  the  individual  instantaneous  values,  V  and  T,  in 
terms  of  their  means  and  anomalies,  as 

V  =  V  +  V'  ,  V  =  [V]  +  V* ,  T  =  T  +  T'  ,  T  =  [T]  +  T* 

The  bars  and  brackets  have  their  previously  defined  meanings,  and  the 
prime  and  asterisk  superscripts  represent  departures  from  the  time  and 
space  averages,  respectively.  By  making  the  above  substitutions,  we  see 
that 

[VT]  =  [V] [Tj  +  [ V*T* ]  +  [V'T']  (1.3) 

recognizing  that  the  time  and  zonal  averages  of  the  departures  are  zero 
by  definition. 

The  first  term  on  the  right  hand  side  of  equation  1.3  represents 
the  heat  transport  due  to  the  mean  meridional  circulation.  This  is  the 
heat  transport  due  to  the  Hadley,  Ferrel,  and  Polar  cells.  The  last  two 
terms  represent  the  meridional  eddy  heat  transport  broken  down  into  its 
components.  The  term  [V*T*]  depicts  the  standing  or  stationary  eddy  heat 
transport.  Stationary  eddies  are  waves  in  the  mean  flow  that  persist  over 
the  averaging  time  period.  For  example,  if  we  select  an  averaging  time  of 
a  month,  then  the  waves  on  a  monthly  mean  500  mb  map  would  be  stationary 
eddies.  The  last  term  in  equation  1.3  represents  the  meridional  transient 
eddy  heat  transports.  Transient  eddies  are  all  deviations  in  the  eddy 
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circulations-  within  the  averaging  time  period.  For  example,  according  to 
Clapp  (1970) ,  if  the  averaging  time  interval  is  a  month,  transient  eddies 
are  the  part  of  the  mean  meridional  heat  transport  due  largely  to  travel¬ 
ing  cyclones  and  anticyclones. 

Steady  state  models  used  to  simulate  climate  cannot  calculate 
transient  eddy  heat  transports.  The  other  two  terms  (the  stationary  eddy 
and  mean  meridional  circulation  heat  transports)  involve  the  covariance 
or  product  of  the  time-averaged  quantities  and  can,  therefore,  (at  least 
in  principle)  be  predicted  explicitly  by  a  steady  state  model  (Clapp,  1970) 
For  this  reason,  studies  have  been  directed  toward  the  understanding  and 
parameterization  of  the  transient  eddy  heat  transports. 

1.2  Brief  Review  of  Previous  Work 

Parameterization  of  transient  eddy  heat  transport  using  an  "Aust- 
ausch  coefficient"  approach  was  originally  proposed  by  Defant  (1921) .  In 
this  approach,  the  transient  eddy  heat  transport  term  in  equation  1.3  is 
approximated  by 

[T'V']  =  -K(9[T]/9y)  (1.4) 

where  the  bar  and  brackets  signify  time  and  zonal  averages  as  before,  and 
K  is  the  Austausch  coefficient.  Vrtiite  and  Jung  (1951)  were  among  the 
first  to  estimate  the  Austausch  coefficient  using  mean  heat  transports 
computed  from  synoptic  weather  maps.  In  their  study,  they  also  observed 
a  negative  correlation  between  eddy  heat  transports  and  temperature 
gradients  for  short  averaging  periods  up  to  twelve  days;  however,  they  did 


s 
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not  separate  stationary  and  transient  eddy  components  in  their  computa¬ 
tions.  Later  Saltzman  (1967) ,  basing  his  ideas  on  the  linear-pertur¬ 
bation  form  of  the  hydrodynamic  and  thermodynamic  equations,  suggested 
that 


K  =  -B(3[T]/3y)  (1.5) 

* 

where  B  is  a  stability  coefficient.  Substituting  equation  1.5  into  equa¬ 
tion  1.4  gives 

('FVM  =  B  (3  [T]/3Y) 2  (1.6) 

indicating  that  the  transient  eddy  heat  transport  is  proportional  to  the 
square  of  the  temperature  gradient. 

Clapp  (1970)  used  two  independent  sets  of  data  to  obtain  new  esti¬ 
mates  of  K  and  B,  and  investigated  the  Austausch  formulae.  His  investiga¬ 
tion  suggests  that  the  Austausch  formulae  may  be  fairly  successful  in  es¬ 
timating  the  zonally  averaged  meridional  transport.  In  attempting  to  ex¬ 
tend  this  test  to  explain  longitudinal  variations,  his  preliminary  efforts 
were  not  successful.  This  suggests  that  further  investigations  and  obser¬ 
vations  of  meridional  eddy  heat  transports  are  needed,  especially  studies 
containing  longitudinal  variations . 

Oort  and  Rasmusson  (1971)  have  compiled  a  very  complete  set  of 
computed  transports  for  the  five  year  period  May  1958  to  April  1963.  Their 
computations  included  meridional  heat  transports  for  each  month  broken  down 
into  the  three  components  described  previously.  All  their  work  was  direct¬ 
ed  at  zonally  averaged  quantities ,  so  they  could  not  observe  longitudinal 


differences.  Latitudinal,  seasonal  and  monthly  variations  in  the  heat 


13  - 


transports  were  observed  and  discussed.  Many  of  their  observations  are 
relevant  to  our  investigation,  and  we  will  refer  to  their  work  often  in 
the  course  of  this  paper. 

Blackmon  et  al.  (1977)  performed  an  observational  and  statistical 
study  of  heat  fluxes.  They  used  time-filters  to  study  fluctuations  of 
different  periods.  The  most  definitive  results  involved  "band-pass"  flue- 
tuations  (fluctuations  having  a  period  of  2.5-6  days)  which  appeared  to  be 
associated  with  developing  baroclinic  waves.  Their  investigation  suggested 
a  relationship  between  poleward  eddy  heat  fluxes  and  baroclinic  instability 
associated  with  the  strong  thermal  gradients  at  the  earth's  surface.  This 
is  one  of  the  results  that  influenced  the  structure  of  our  study. 

The  main  motivation  for  our  study  was  a  paper  by  Stone  (1978) .  In 
his  paper,  he  compared  zonal  mean  meridional  temperature  gradients  in  the 
atmosphere  to  critical  temperature  gradients  predicted  by  a  two- layer  baro¬ 
clinic  model.  Stone  observed  that  eddy  heat  fluxes  are  sensitive  to  changes 
in  meridional  temperature  gradients  (i.e.  baroclinic  instability).  This 
observation  suggests  a  relationship  between  these  two  quantities. 

We  intend  to  statistically  investigate  the  relationship  between 
meridional  eddy  heat  transports,  meridional  surface  temperature  gradients, 
and  the  two-layer  model  stability  parameter,  since  baroclinic  theory  and 
results  from  previous  studies  have  suggested  a  possible  relationship  be¬ 
tween  these  quantities.  We  will  deal  only  with  the  Northern  Hemisphere, 
where  the  raw  data  are  more  abundant.  In  order  to  observe  latitude  and 
longitude  variations,  quantities  were  used  which  had  been  calculated  for 
evenly  spaced  grid  points.  In  this  study,  values  of  the  total  tropospheric 
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meridional  eddy  heat  transport  were  used  which  were  defined  as  follows: 


EF 


h 

[  PC  (V-[V]) (T-[T])dZ 


(1.7) 


A  time  averaging  period  of  a  month  was  used,  and  h  is  the  height  of  the 
tropopause.  This  flux  was  divided  into  two  components.  The  meridional 
stationary  eddy  heat  transport  at  each  grid  point  was  computed  from  the 
equation 

h 

SE  =  |  pCp(V*T*)dZ  (1.8) 

O 

i.e.,  SE  represents  the  monthly  mean  meridional  stationary  eddy  heat  trans¬ 
port.  Once  the  stationary  and  total  eddy  heat  transport  were  computed,  a 
transient  eddy  heat  transport  was  simply  determined  by  subtraction. 


TE  *  EF  -  SE 


(1.9) 


where  TE  is  the  monthly  mean  meridional  transient  eddy  heat  transport.  It 
is  important  to  note  that  fluctuations  in  the  meridional  mean  circulation 
are  not  included  in  the  TE  term,  but  are  included  in  the  transient  eddy 
heat  transport  term  in  equation  1.3. 

With  the  quantities  calculated  at  each  grid  point,  we  are  able  to 
perform  a  correlation  analysis  for  all  possible  combinations  of  these  quan¬ 
tities.  We  will  concentrate  our  investigation  on  longitudinal  variations, 
since  this  area  has  not  been  explored.  Once  the  correlation  analysis  is 
performed,  empirical  parameterization  schemes  for  the  meridional  transient 
eddy  heat  transport  will  be  investigated. 


CHAPTER  TWO 


SOURCES  AND  DESCRIPTION  OF  DATA 

2.1  Meridional  Eddy  Heat  Transport  Data 

The  meridional  eddy  heat  transport  data  was  calculated  at  the 
Goddard  Institute  for  Space  Studies  (GISS) ,  under  the  direction  of  Pro¬ 
fessor  Peter  H.  Stone,  Massachusetts  Institute  of  Technology.  The  raw 
data  for  the  calculations  was  provided  by  the  National  Meteorological  Cen¬ 
ter  (NMC) .  Mean  monthly  values  of  eddy  heat  transport  were  available  for 
January  1973,  1974,  and  1975.  The  data  was  recorded  for  every  four  degrees 
of  latitude  90°S,  86°S,  ...  86°N,  90°N,  and  every  ten  degrees  of  longitude 
175°W,  165°W,  ...  165°E,  175°E.  This  study  was  based  on  the  data  between 

18°N  and  70°N  latitude.  All  the  meridional  eddy  heat  transport  data  are 

17 

vertically  averaged  values  for  the  troposphere,  and  are  in  units  of  10 
calories  per  day  per  five  degrees  longitude  with  positive  indicating 
northward.  Before  the  vertical  averaging,  the  tropopause  was  determined 
separately  for  each  grid  point.  The  meridional  eddy  heat  transport  data 
was  separated  into  stationary  eddy  and  transient  eddy  components  for  the 
monthly  mean  data. 

2.1.1  Meridional  Stationary  Eddy  Heat  Transport  Data 
Oort  and  Rasmusson  (1971)  observed  that  the  meridional  stationary 
eddy  heat  transports  (MSEHT)  were  strongest  in  the  winter.  This  is  one 
of  the  main  reasons  the  month  of  January  was  used  for  our  study.  Haines 
and  Winston  (1963)  first  noted  that  the  MSEHTs  were  dominated  by  three 
main  features,  the  Aleutian  low,  the  Icelandic  low,  and  the  Siberian  high 


-  16  - 


pressure  systems.  They  also  observed  that  the  peaks  of  the  transports  are 
located  to  the  east  of  these  features,  i.e.  the  transport  of  cold  air  south¬ 
ward  by  the  northerly  flow  over  eastern  Siberia  and  the  transport  of  warm 
air  northward  over  the  central  North  Atlantic  and  Gulf  of  Alaska.  The 
exact  position  and  relative  importance  of  these  features  varies  from  year 
to  year,  as  pointed  out  by  Blackmon  et  al.  (1977).  Tables  2.1,  2.2,  and 
2.3  contain  the  MSEHT  data  for  our  three  months.  From  the  tables,  we  can 
see  that  the  Icelandic  low  was  very  strong  and  the  dominant  feature  for 
the  monthly  mean  MSEHT  in  January  1973.  For  January  1974,  all  three  fea¬ 
tures  were  strong,  but  the  Icelandic  low  was  still  slightly  dominant  and 
shifted  20°E  of  its  previous  years  position.  In  January  1975,  the  Siberian 
high  and  the  Aleutian  low  pressure  systems  were  the  dominant  features. 

The  latitude  of  the  peak  MSEHT  fluctuates  annually.  Table  2.4  con¬ 
tains  the  MSEHT  summed  around  each  latitude  circle.  By  examining  the  summed 
values  of  our  MSEHT,  we  see  the  peak  transport  occurred  at  58°N  for  January 
1973.  In  January  1974,  the  peak  MSEHT  was  located  at  46°N.  January  1975 's 
peak  was  situated  at  50°N.  These  fluctuations  from  year  to  year  indicate 
not  only  the  variability  of  the  transports,  but  also  that  we  have  a  variety 
of  data  to  study,  rather  than  a  bias  sampling  of  one  particular  case.  The 
position  of  the  peak  MSEHT  for  the  average  over  the  three  months  was  at 
latitude  50°N.  Oort  and  Rasmusson  (1971)  also  found  the  peak  MSEHT  in  this 
location  with  data  over  a  five  year  period. 

2.1.2  Meridional  Trams ient  Eddy  Heat  Tramsport  Data 


The  meridional  transient  eddy  heat  transport  (MTEHT)  is  of  the  same 
order  of  magnitude  as  the  MSEHT,  when  comparing  the  summed  values  in  Table 
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2.4;  however,  the  peaks  of  the  MSEHT  are  generally  larger  than  the  peaks  of 
the  MTEHT.  The  MTEHT  are  spread  more  uniformly  over  the  globe  than  the 
MSEHT.  Even  though  the  MTEHTs  are  dispersed  around  the  globe, 
there  are  locations  where  peak  transports  do  exist.  Tables  2.5, 
2.6,  and  2.7  display  the  transient  compound  of  the  meridional 
eddy  heat  transports.  Notice  that  the  position  of  the  peaks 
supports  the  observation  made  by  Blackmon  et  al .  (1977)  that  the 

peaks  are  closely  related  to  the  major  storm  tracks.  Three  major 
storm  tracks,  along  the  east  coast  of  the  United  States  extending 
to  Greenland,  along  the  east  coast  of  Asia  up  to  Alaska,  and  a 
short  track  along  the  west  coast  of  the  United  States  and  Canada , 
are  primarily  emphasized  by  our  data.  The  exact  position  and 
strength  of  these  peaks  varies  annually.  In  January  1973,  the 
peaks  were  mainly  along  the  east  coasts  of  Asia  and  the  United 
States.  These  two  peaks  were  again  the  dominant  features  in 
January  1974;  however,  they  were  displaced  further  north  and  east 
of  the  previous  year's  position.  January  1975  was  dominated  by 
three  peaks.  The  peaks  along  the  east  coasts  of  Asia  and  the 
United  States  were  still  evident  and  were  located  primarily  be¬ 
tween  the  January  1973  and  January  1974  positions.  A  third  peak 
was  located  along  the  west  coast  of  Canada  and  the  United  States. 
It  was  smaller  in  area  coverage  but  of  the  same  magnitude  as  the 
other  two  peaks.  The  variability  of  these  peaks  reflect  the 
variability  of  the  storm  tracks  from  year  to  year. 

Comparing  the  summed  values  of  the  MSEHT  and  the  MTEHT  listed  in 
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Table  2.4,  we  can  see  the  annual  variability  of  the  dominance  of  the  sta¬ 
tionary  and  transient  components  with  latitude.  The  transient  component 
was  larger  at  all  the  latitudes  for  January  1975.  Just  the  opposite  was 
observed  in  January  1974.  January  1973  had  the  stationary  component  larger 
from  50°N  to  70°N,  and  the  transient  component  was  larger  from  30°N  to  46°N. 

Also  looking  at  Table  2.4,  we  see  that  the  peak  of  the  MTEHT  varied 
in  latitude  from  year  to  year.  In  January  1973,  the  peak  was  at  42#N. 

The  peak  in  January  1974  was  at  50°N.  Latitude  46°N  was  the  location  of 
the  peak  in  January  1975.  When  averaged  over  the  three  months,  the  peak 
was  at  46°N.  Oort  and  Rasmusson  (1971)  also  observed  the  peak  in  this  lo¬ 
cation  over  their  five  years  of  data.  This  indicates  our  data  set  is  rep¬ 
resentative  of  a  typical  set  of  Januaries. 

2.1.3  Total  Meridional  Eddy  Heat  Transport  Data 

Since  the  total  meridional  eddy  heat  transport  is  just  the  sum  of 
the  stationary  and  transient  components,  it  contains  some  of  the  charac¬ 
teristics  described  above.  These  transports  are  displayed  in  Tables  2.8, 
2.9,  and  2.10.  The  total  meridional  eddy  heat  transport  resembles  the 
MSEHT  to  a  large  extent,  because  of  the  dominance  of  the  stationary  peaks. 
This  is  particularly  true  because  there  are  some  similarities  between  the 
stationary  and  transient  components  so  that  when  added  together  they  rein¬ 
force  each  other.  This  is  especially  noticeable  in  the  North  Atlantic. 

Looking  at  the  summed  values  in  Table  2.4,  we  see  that  the  peak 
transports  did  not  vary  with  latitude  annually.  For  all  three  months, 
the  peak  was  found  at  50 °N.  Oort  and  Rasmusson  (1971)  observed  that  over 
their  five-year  period,  the  peak  transport  was  also  located  at  50°N. 


Table 
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2.2  Meridional  Surface  Temperature  Gradient 

2.2.1  Acquisition  of  Data 

Land  stations'  mean  surface  temperatures  for  January  1973,  1974, 
and  1975  were  obtained  from  "Climatic  Data  of  the  World",  published  month¬ 
ly  by  the  U.S.  Environmental  Data  Service  (NOAA) .  Effects  of  elevation 
were  taken  into  account  by  reducing  the  temperature  to  sea  level  based  on 
the  U.S.  Standard  Atmosphere  lower  tropospheric  lapse  rate  of  6.5  degrees 
Celsius  per  kilometer.  The  temperatures  were  plotted  and  analyzed  on  a 
northern  hemispheric  map.  Temperatures  every  ten  degrees  of  longitude  and 
every  four  degrees  of  latitude  were  read  and  recorded  on  a  separate  table 
to  use  m  obtaining  meridional  temperature  gradients  at  the  points  where 
the  meridional  eddy  heat  transport  data  was  available.  Meridional  surface 
temperature  gradients  were  determined  by  taking  north-south  centered  dif¬ 
ferences,  and  are  recorded  in  units  of  degrees  Celsius  per  four  degrees 
latitude.  Negative  values  mean  that  the  temperature  is  decreasing  toward 
the  north.  The  meridional  surface  temperature  gradients  are  listed  in 
Tables  2.11,  2.12,  and  2.13. 

Some  problems  with  the  data  occurred,  causing  the  analysis  to  be 
somewhat  subjective.  High  elevation  stations,  when  reduced  to  sea  level, 
gave  unrepresentative  temperatures  compared  to  the  surrounding  lower  ele¬ 
vation  stations.  The  problem  was  primarily  in  the  Alps,  and  these  stations 
were  discarded.  Large  data-sparse  areas  occurred  in  the  Sahara,  the  Arab¬ 
ian  Peninsula,  and  the  People's  Republic  of  China.  In  these  areas,  we 
followed  the  general  pattern  of  the  surrounding  isotherms  and  thus  lost 
any  small  scale  features  which  may  have  existed. 
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Air  temperature  data  was  not  available  over  the  oceans,  and  so  we 
used  sea  surface  temperatures.  Although  these  temperatures  may  not  relate 
directly  to  corresponding  air  temperatures,  their  gradients  are  likely  to 
correspond.  Mean  sea  surface  temperatures  were  available  directly  from 
data  compiled  by  the  National  Marine  Fisheries  Service  (NOAA) .  The  tem¬ 
peratures  were  averaged  values  over  a  five-by-five  degree  square  and  were 
in  degrees  Celsius.  They  were  plotted,  analyzed,  and  the  gradients  were 
determined  in  the  same  manner  as  over  the  land.  Whenever  land  and  sea  data 
analysis  overlapped,  the  land  analysis  was  used.  This  did  not  occur  enough 
to  make  a  reasonable  comparison  of  the  gradients. 

2.2.2  Data  Characteristics 

The  strongest  meridional  surface  temperature  gradients  for  all 
three  months  occurred  where  oceans  lie  to  the  south  of  land  in  the  upper 
mid-latitudes.  The  strongest  of  these  was  in  eastern  Asia,  north  of  Japan. 
Southern  Alaska  and  western  Canada's  coastline  also  had  a  very  strong  sur¬ 
face  temperature  gradient  that  was  present  for  all  three  months.  Looking 
at  the  zonally  averaged  meridional  surface  temperature  gradients,  we  see 
that  the  peak  was  always  at  62°N.  In  January  1974  and  1975,  the  tempera¬ 
ture  gradient  on  the  south  coast  of  Nova  Scotia  was  quite  noticeably  strong; 
however,  in  January  1973  the  temperature  gradient  was  only  half  as  intense. 
Minor  peak  temperature  gradients  occurred  just  to  the  south  of  the  semi¬ 
permanent  cold  core  highs  in  Siberia  and  Northern  Canada.  The  Siberian 
high  produced  the  stronger  of  the  two  gradients.  The  position  and  magni¬ 
tude  of  these  last  two  meridional  surface  temperature  gradients  varied 
slightly  from  year  to  year.  The  data  south  of  30°N  should  not  be  trusted 
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to  any  large  extent,  because  this  area  contains  large  data-sparse  areas, 
and  very  large  areas  where  sea  surface  temperature  data  was  used. 

2.3  Two-Layer  Model  Stability  Parameter 

In  the  two-layer  model,  the  troposphere  is  divided  into  two  equal 

mass  layers.  The  mass  averaged  parameters  describing  the  lower  layer  will 

be  denoted  by  the  subscript  1,  the  subscript  2  will  denote  the  upper  layer. 

The  stability  parameter  can  now  be  defined  as:  U_  -  U,  -  U  j  Stone  (1978). 

U_  -  U,  is  just  the  vertical  shear  of  the  two  layer  troposphere.  U  is 
2  1  c 

the  critical  wind  shear,  and  was  determined  by  Phillips  (1954)  and  Stone 
(1978)  a  . 


U 


BR(02-ei) 


(2.1) 


where  B  =  ,  f  =  2  sin  <t>,  R  is  the  gas  constant,  and  0j  and  02  are  the 

mass  averaged  potential  temperatures  of  the  two  layers.  The  stability 
parameter  was  provided  from  the  same  source  and  in  the  same  format  as  the 
meridional  eddy  heat  transport  data.  Units  are  meters  per  second,  and 
positive  values  indicate  instability,  while  negative  values  indicate  sta¬ 
bility. 


2.3.1  Model  Stability  Parameter  Characteristics 

The  model  stability  parameter  is  noisy,  mainly  because  of  the  large 
variation  in  U2  -  U^.  To  solve  this  problem,  the  stability  parameter  was 
smoothed  over  a  number  of  degrees  of  longitude.  The  smoothing  distance 
was  2Lr<$)  rounded  to  the  nearest  5  degrees,  where  Lr ($)  is  the  mean  radius 
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of  deformation.  For  example,  the  smoothing  distance  for  the  latitudes  38°N 
to  SO°N  was  40  degrees  longitude. 

To  a  good  approximation,  was  constant  along  a  latitude  circle 
(Stone,  personal  communication) .  Having  U essentially  constant  on  a  lati¬ 
tude  circle  enables  us  to  calculate  the  variation  of  the  vertical  wind 
shear  from  the  variation  of  the  model  stability  parameter.  The  vertical 
wind  shear  can  be  related  to  the  horizontal  temperature  gradient  by  the 
thermal  wind  relationship.  Thus  to  a  large  extent,  we  will  be  correlating 
the  mean  meridional  tropospheric  horizontal  temperature  gradient  to  the 
meridional  eddy  heat  transports,  when  we  correlate  the  latter  with  the 
stability  parameter. 

2.3.2  Data  Characteristics 

Tables  2.14,  2.15,  and  2.16  contain  the  stability  parameter.  One 

striking  feature  of  the  data  was  the  dominance  of  negative  values  south  of 

30°N.  The  reason  for  this  is  two-fold.  Primarily,  it  is  because  U  increases 

c 

as  the  latitude  decreases.  Also,  the  vertical  wind  shear  is  much  smaller 
in  the  tropics  than  in  the  extratropical  latitudes,  so  decreases  as 

we  move  into  the  tropics.  Obviously,  the  atmosphere  is  not  as  completely 
stable  as  the  model  stability  parameter  indicates  at  and  below  30°N,  since 
the  eddy  transports  are  not  zero.  Nevertheless,  they  are  relatively  small 
in  this  region. 

There  are  two  main  areas  where  the  mean  monthly  stability  parameter 
has  a  large  positive  value  for  all  three  months.  One  unstable  area  is  just 
off  the  east  coast  of  Japan.  The  other  unstable  area  is  located  south  of 
Nova  Scotia.  It  is  interesting  to  notice  that  these  areas  are  located  in 
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the  vicinity  of  the  major  winter  storm  tracks.  The  position  and  magnitude 
of  these  unstable  region  varies  slightly  from  year  to  year.  The  area  off 
the  east  coast  of  Japan  had  the  largest  values  for  January  1973,  and  1975. 
The  area  south  of  Nova  Scotia  was  the  dominant  peak  in  January  1974.  For 
January  1973,  both  areas  were  noticeably  weaker  compared  with  the  other 
two  months.  A  third,  smaller  unstable  area  appeared  in  January  1975  along 
the  west  coast  of  Canada.  This  area  was  not  present  in  the  previous  two 
years. 
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CHAPTER  THREE 
METHODS  OF  ANALYSIS 


3.1  Correlation  Analysis 

Correlations  were  only  performed  along  the  latitude  circles  for 
several  reasons.  Latitudinal  correlations  were  not  performed  because  of 
the  limited  degrees  of  freedom  in  that  direction,  and  because  the  results 
are  fairly  obvious  just  by  examining  the  data  tables  (i.e.  peak  transports 
occur  in  mid-latitudes  where  baroclinic  instability  is  the  largest) .  Also, 
latitudinal  correlations  were  not  computed  because  the  quality  of  the  data 
was  not  as  good  below  34°N  latitude.  Therefore,  each  of  the  fourte''  .  lat¬ 
itudes  were  treated  separately.  It  is  important  to  note  that  any  similar¬ 
ities  between  different  latitude  correlations  must  be  attributed  solely  to 
the  data,  so  the  persistence  of  large  correlations  at  different  latitudes 
adds  to  the  significance  of  the  correlations. 

The  correlation  coefficients  for  the  36  points  around  the  latitude 
circle  was  computed  as  follows: 

r  =  E  (3.1) 

36  0X  Oy 

where  r  is  the  correlation  coefficient,  X  and  Y  are  the  variable  valves, 

Xgj  and  Yffl  are  the  zonal  averages  of  the  two  variables,  and  0X  and  CTy  are 
the  individual  standard  deviations.  The  correlation  coefficients  were  not 
only  computed  directly,  but  also  using  spatial  lags  of  10  to  350  degrees 
of  longitude  in  increments  of  ten  degrees.  To  better  explain  the  use  of 


* 
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the  term  spatial  lags,  let  us  take  an  example  of  correlating  two  arrays  X 
versus  Y.  This  means  that  we  are  correlating  the  first  array  X(A)  with 
the  second  array  Y (A+A) ,  where  A  is  the  longitude,  and  A  =  0,  10,  20,  ... 
350  degrees  longitude  is  the  eastward  shift  of  the  second  variable  (spa¬ 
tial  lag) .  This  allowed  us  to  calculate  36  correlation  coefficients  for 
each  of  the  14  latitudes,  for  each  pair  of  input  arrays.  Since  the  data 
along  the  latitude  circle  looped  completely  around  the  earth,  there  was  no 
lop-off  error  with  the  spatial  lags. 

Each  January  had  five  input  arrays.  By  correlating  each  input 
array  with  itself  and  the  other  four  arrays,  we  obtained  fifteen  different 
correlation  coefficient  arrays  per  year.  Of  the  fifteen  combinations,  five 
were  autocorrelations  which  were  used  for  determining  the  degrees  of  free¬ 
dom  in  the  data.  The  other  ten  combinations  are  the  crux  of  our  investi¬ 
gation.  By  comparing  these  ten  different  combinations  for  each  January, 
it  was  obvious  that  the  patterns  were  quite  similar.  There  were  some  minor 
differences  from  year  to  year,  and  we  will  refer  to  these  differences  as 
noise.  This  noise  may  have  been  caused  by  any  of  a  number  of  factors, 
such  as  errors  in  data  measurements  or  handling,  or  natural  variability, 
etc.  Our  data  does  not  span  a  long  enough  period  to  test  any  explanation 
for  the  interannual  variability,  so  the  noise  is  of  little  importance  to 
this  investigation.  To  reduce  the  noise  level,  we  averaged  the  correla¬ 
tion  coefficient  arrays  over  the  three  Januaries.  The  averaging  reduced 
the  peak  correlation  coefficients,  but  did  not  change  the  significance  of 
the  important  peaks,  because  it  tripled  the  number  of  data  points,  thus 
increasing  the  degrees  of  freedom.  These  averaged  correlation  coefficient 
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arrays  are  contained  in  the  tables  and  will  be  discussed  in  the  next  chapter. 


3.2  Significance  Testing 

The  major  difficulty  in  significance  testing  is  determining  the 
degrees  of  freedom  of  the  data.  To  determine  the  degrees  of  freedom,  we 
used  a  method  prescribed  by  Davis  (1976) .  In  our  case,  the  degrees  of 
freedom  (N)  are  defined  as 


N 


where 


(3.2) 


An 

AA 


(3.3) 


C  (A{ )  and  C  (A^ )  are  the  autocorrelation  functions  for  the  two  correlated 
x  1  y  1 

data  arrays  averaged  over  the  three  Januaries.  The  three  in  the  definition 
of  N  is  present  because  we  are  averaging  over  the  three  completely  inde¬ 
pendent  Januaries.  Unlike  Davis,  we  limited  our  spatial  lags  to  +  90  de¬ 
grees  in  order  to  avoid  the  influence  of  continents,  which  introduces  a 
strong  wave  number  two  component  in  the  data.  AS  it  turned  out,  this  pro¬ 
cedure  only  slightly  affected  our  confidence  level  determination.  The 
An 

series  for  converges  very  rapidly  as  we  can  see  in  the  examples  in 

Figure  3.1.  Because  of  the  continental  effect,  we  limit  our  discussion  of 
correlation  peaks  to  those  within  ±  90  degrees  spatial  lag. 

Once  the  degrees  of  freedom  were  known,  it  was  a  simple  procedure 


to  use  the  t  test  tables  to  determine  the  95%  and  99%  confidence  levels 
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for  an  individual  correlation.  These  confidence  levels  apply  to  a  priori 
correlation  peaks.  A  priori  correlation  peaks  are  correlation  peaks  cen¬ 
tered  around  zero  degrees  spatial  lag,  because  correlation  peaks  at  zero 
degrees  spatial  lag  are  physically  expected.  Correlation  peaks  centered 
at  other  spatial  lags  are  a  posteriori  peaks.  We  are  limited  in  our  abil¬ 
ity  to  test  the  significance  of  the  a  posteriori  peaks,  because  we  only 
have  three  months  of  data.  This  allows  us  only  three  degrees  of  freedom 
for  these  peaks,  giving  their  95%  and  99%  confidence  levels  as  .80  and  .93, 
respectively.  A  larger  data  base  would  be  necessary  to  fully  test  the 
relationship  suggested  by  peaks  at  finite  lags. 

The  confidence  levels  for  the  a  priori  correlation  peaks  ar  listed 
in  Table  3.1  for  the  two  correlations  of  most  interest.  There  was  no  sig¬ 
nificance  testing  performed  for  the  other  combinations,  since  preliminary 
investigation  showed  results  of  little  importance  (Stailey,  personal  com¬ 
munication).  It  is  interesting  to  note  from  Table  3.1  that  there  is  little 
latitudinal  variation  in  the  99%  confidence  level  for  the  combination  of 
surface  temperature  gradient  versus  transient  eddy  heat  transport.  Lati¬ 
tudinal  variation  is  a  little  more  pronounced  for  the  model  stability  para¬ 
meter  versus  transient  eddy  heat  transport,  and  the  degrees  of  freedom  are 
noticeably  smaller  than  for  the  first  combination.  This  can  be  explained 
by  the  fact  that  the  model  stability  parameter  was  smoothed  longitudinally. 
One  final  point  to  note  is  that  the  degrees  of  freedom  are  fewest  in  mid¬ 
latitudes  for  both  combinations.  This  would  indicate  that  the  variables 
do  not  differ  as  rapidly  with  longitude  as  they  do  for  the  other  latitude 


% 


circles . 
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CHAPTER  FOUR 

DISCUSSION  OF  CORRELATIONS 

All  of  the  averaged  correlation  coefficient  arrays,  except  for  the 
autocorrelations  and  the  two  trivial  combinations,  transient  eddy  heat 
transport  versus  total  eddy  heat  transport  and  stationary  eddy  heat  trans¬ 
port  versus  total  eddy  heat  transport,  will  be  discussed  in  this  chapter. 
The  discussion  will  be  confined  primarily  to  the  area  between  and  including 
latitudes  34°N  to  70°N  and  +  90  degrees  spatial  lag,  for  the  reasons  men¬ 
tioned  earlier.  Significance  testing  was  performed  on  the  two  combina¬ 
tions,  surface  temperature  gradient  vr  sus  transient  eddy  heat  transport 
and  two-layer  model  stability  parameter  versus  transient  eddy  heat  trans¬ 
port,  because  from  baroclinic  theory  one  expects  a  correlation.  For  these 
combinations,  the  correlation  coefficients  greater  than  or  equal  to  the 
99%  confidence  level  for  an  individual  correlation  are  contoured  with  a 
solid  line.  The  other  combinations  are  also  of  some  interest,  and  for 
these  combinations,  correlation  coefficients  exceeding  .34  are  contoured 
with  a  dashed  line.  The  .34  cutoff  was  selected  because  this  would  ap¬ 
proximately  be  the  99%  confidence  level. 

4.1  Meridional  Surface  Temperature  Gradient  versus  Two-Layer  Model 
Stability  Parameter 

The  results  are  contained  in  Tables  4.1A  and  4.1B.  The  large 
negative  correlations  are  easily  explained.  Since  Uc  is  essentially  in¬ 
dependent  of  longitude,  it  is  normalized  out  in  the  correlation  process, 
and  in  essence  we  are  correlating  the  meridional  surface  temperature  gra¬ 
dient  with  the  meridional  mean  tropospheric  temperature  gradient.  It  is  not 


SURFACE  TEMPERATURE  GRADIENT  VS  MODEL  STABILITY  PARAMETER  (JAN  73#  74»  75) 


Table  ij.,1  A 
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surprising,  then,  to  have  a  strong  correlation  between  these  two  quantities  near 
zero  degrees  spatial  lag.  The  negative  sign  exists  because  of  the  opposite  sign 
convention  of  the  two  quantities .  A  large  negative  surface  temperature  gradient 
means  a  strong  northward  decrease  in  temperature;  whereas  a  large  positive  model 
stability  parameter  means  large  baroclinic  instability,  which  in  the  Northern 
Hemisphere  indicates  a  strong  decrease  in  temperature  to  the  North. 

The  strongest  negative  correlations  occurred,  on  the  average,  at  a 
spatial  lag  of  10  to  20  degrees.  This  indicates  that  the  peak  baroclinic 
instability  occurs  10  to  20  degrees  downstream  (to  the  east)  of  the  peak 
surface  temperature  gradients.  Obviously,  this  just  reveals  the  vertical 
slope  of  the  mean  temperature  field.  Tt  is  interesting  to  observe  that 
the  only  latitude  which  strongly  deviates  from  this  behavior  is  latitude 
62°N,  where  the  surface  temperature  gradients  were  the  strongest. 

Looking  at  the  spatial  lags  270  to  330  degrees,  we  see  positive 
peak  correlations.  These  correlations  appear  to  differ  in  location  with 
latitude;  however,  on  the  average  the  positive  peak  correlations  are  about 
90  degrees  apart  from  the  negative  peaks.  This  indicates  that  both  quan¬ 
tities  have,  in  general,  a  two  wave  cycle  around  the  earth,  corresponding 
with  two  major  land  masses  and  two  oceans.  The  variation  in  latitude  can 
be  explained  by  the  variation  in  land  and  ocean  distribution  around  the 
latitude  circles.  At  and  above  58  °N,  the  latitude  circles  are  dominated 
by  land  masses.  This  would  explain  the  disruption  in  the  two  wave  pattern 
at  these  latitudes.  The  two  wave  pattern  is  clearly  visible  in  the  lati¬ 
tudes  from  42°N  to  54°N.  South  of  42°N,  the  continent  of  North  America 
becomes  much  smaller  than  the  Asian  continent,  and  again  the  two  wave  pat¬ 
tern  is  disrupted.  Looking  at  the  autocorrelations,  this  continental 


effect  is  most  strongly  evident  with  the  surface  temperature  gradients. 

Interannual  variability  is  interesting  in  this  case.  The  corre¬ 
lation  patterns  and  values  are  very  similar  for  January  1973  and  1974. 

The  peak  correlation  coefficients  occurred  at  latitude  58°N,  with  spatial 
lags  of  10  and  20  degrees.  The  average  correlation  coefficient  for  this 
location  was  -.74.  January  1975  was  radically  different,  especially  for 
latitude  58°N.  The  correlations  were  weaker  and  shifted  considerably 
farther  eastward.  To  explain  this  phenomenon,  we  would  need  a  much  larger 
data  base,  but  this  does  suggest  that  the  relationship  between  these  two 
quantities  is  dependent  on  other  variables . 

4.2  Meridional  Stationary  Eddy  Heat  Transport  versus  Meridional  Transient 
Eddy  Heat  Transport 

The  relationship  between  the  two  components  of  the  total  meri¬ 
dional  eddy  heat  transport  was  investigated  by  using  the  correlation  analy¬ 
sis  scheme.  Tables  4.2A  and  4.2B  contain  the  results.  Looking  at  these 
tables,  we  see  there  is  no  strong  systematic  relationship  between  these 
two  components.  We  do  have  a  peak  in  the  correlation  coefficients  at 
latitudes  34°N  and  38°N  with  spatial  lags  of  10  and  20  degrees,  but  this 
peak  is  weak  and  covering  a  very  few  points,  suggesting  it  is  not  signifi¬ 
cant.  Another  peak  is  centered  at  latitude  62°N  with  a  spatial  lag  of  200 
degrees.  The  spatial  lag  is  so  large  that  the  significance  of  this  peak 
is  questionable  at  best.  Therefore,  no  obvious  relationship  exists  between 
these  components.  Looking  at  the  individual  monthly  correlations,  we  see 
this  observation  supported  in  each  case.  This  suggests  that  the  position 
and  strength  of  the  stationary  and  transient  eddy  heat  transports  are 
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mainly  caused  by  different,  and  unrelated  factors. 

4.3  Meridional  Surface  Temperature  Gradient  versus  Meridional  Stationary 
Eddy  Heat  Transport 

Due  to  the  sign  convention  used  for  the  surface  temperature  gra¬ 
dient,  a  direct  relationship  will  be  indicated  by  a  negative  correlation. 
This  fact  should  be  remembered  whenever  we  are  discussing  a  correlation 
with  the  surface  temperature  gradient.  Obviously,  Tables  4.3a  and  4.3B 
indicate  that  there  is  no  convincing  relationship  between  the  surface 
temperature  gradient  and  the  stationary  eddy  heat  transport.  We  do  have 
two  peaks  in  our  correlation  coefficients;  however,  they  involve  very  few 
points  and  are  weak.  The  low  correlation  coefficients  cannot  be  attri¬ 
buted  to  the  averaging  because  the  individual  monthly  correlation  coeffi¬ 
cients  are  small  and  there  is  no  large  interannual  variation  in  the  cor¬ 
relation  pattern.  In  general,  it  is  interesting  to  observe  that  the  cor¬ 
relation  coefficients  at  zero  degrees  spatial  lag  are  very  weak.  This 
indicates  that  for  the  monthly  mean  values  the  surface  temperature  gra¬ 
dient  peaks  do  not  occur  in  the  same  location  as  the  stationary  eddy  heat 
transport  peaks.  A  possible  explanation  for  this  observation  is  that  the 
stationary  eddy  heat  transports  are  so  strong  and  effective  in  alleviating 
temperature  differences  that  strong  surface  temperature  gradients  cannot 
develop  in  their  location. 

4.4  Two-Layer  Model  Stability  Parameter  versus  Meridional  Stationary 
Eddy  Heat  Transport 

Tables  4.4A  and  4.4B  contain  the  correlation  analysis.  There  is 
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one  peak  in  the  correlation  coefficients  at  latitude  38 °N  and  a  spatial  lag 
of  20  degrees.  This  peak  contains  only  five  points  and  the  correlation 
coefficients  are  weak,  indicating  that  the  importance  of  the  peak  is  ques¬ 
tionable.  Looking  at  the  individual  monthly  correlation  arrays,  we  noticed 
that  this  peak  was  always  present.  Other  them  this  peak,  no  large  correla¬ 
tions  appeared  in  the  analysis.  It  is  interesting  to  observe  that  this  peak 
is  located  in  the  same  vicinity  as  the  peak  in  the  correlation  analysis  for 
the  stationary  eddy  heat  transport  versus  transient  eddy  heat  transport. 
Table  4.2A.  This  indicates  that  the  only  place  where  the  model  stability 
parameter  and  the  stationary  eddy  heat  transport  is  even  slightly  related  is 
where  the  '•tationary  eddy  heat  transport  behaves  similarly  to  the  transient 
eddy  heat  transport.  Since  the  peak  in  the  correlation  analysis  covers  such 
a  small  number  of  points  and  the  correlation  coefficients  are  weak,  no  strong 
relationship  is  indicated  between  these  two  quantities. 

4.5  Meridional  Surface  Temperature  Gradient  versus  Meridional  Transient 
Eddy  Heat  Transport 

Since  our  correlation  involves  the  surface  temperature  gradients,  a 
direct  relationship  is  indicated  by  a  negative  correlation  coefficient. 
Tables  4.5a  and  4.5B  display  a  large  area  of  high  correlation  coefficients, 
located  between  latitudes  46°N  and  70°N,  and  centered  around  30  to  50  de¬ 
grees  spatial  lag.  The  largest  correlation  coefficient,  -.73,  is  located 
at  latitude  62 °N  and  50  degrees  spatial  lag.  For  the  individual  monthly 
correlation  arrays,  the  strongest  correlation  was  always  found  in  this  loca¬ 
tion.  The  largest  correlation  coefficient, -.83,  was  found  in  January  1973. 
Interannual  variability  of  the  correlation  coefficients  is  quite  interesting 
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in  this  case.  Although  the  largest  correlation  coefficient  always  occurred 
in  the  same  location,  the  area  of  high  correlation  coefficients  varied  from 
year  to  year.  In  January  1973,  the  area  extended  from  latitude  46°N  to 
66°N.  High  correlation  coefficients  existed  only  at  latitudes  62°N  and 
66°N  in  January  1974.  No  high  correlation  coefficients  existed  at  lati¬ 
tude  66°N,  but  they  did  occur  for  the  latitudes  50°N  to  62°N  in  January 
1975.  Latitude  62°N  was  the  only  latitude  where  high  correlation  coeffi¬ 
cients  existed  for  all  three  months.  An  interesting  observation  is  that 
the  strongest  correlation  always  occurred  at  the  latitude  where  the  sur¬ 
face  temperature  gradient  had  its  largest  value.  On  the  other  hand,  where 
the  transient  eddy  heat  transport  was  the  strongest,  at  latitudes  42°N  and 
46°N,  the  correlation  was  weak.  This  would  indicate  that  the  relationship 
between  these  two  variables  is  weak,  and  the  transient  eddy  heat  transports 
may  depend  mainly  on  other  factors.  The  positive  correlation  coefficients 
occurring  in  Table  4.5B  are  due  to  the  wave  characteristics  of  the  two 
quantities  around  the  latitude  circles. 

The  contoured  area  of  correlation  coefficients  in  Table  4.5a  is 
significant  at  the  99%  confidence  level  for  an  a  priori  correlation  peak; 
however,  the  area  exists  away  from  the  zero  degrees  spatial  lag.  This  in¬ 
dicates  that  it  is  am  a  posteriori  probability  peak,  and  has  only  three 
degrees  of  freedom.  In  that  case,  the  95%  confidence  level  is  for  corre¬ 
lation  coefficients  greater  than  .80.  None  of  the  correlation  coefficients 
are  this  large.  The  largest  correlation  is  only  at  the  91%  confidence 
level.  This  does  not  mean  that  the  spatial  relationship  does  not  exist, 
it  just  indicates  that  the  spatial  relationship  is  not  certain.  In  order 
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to  completely  test  this  relationship,  a  larger  data  base  would  be  needed. 

4.6  Two-Layer  Model  Stability  Parameter  versus  Meridional  Transient 
Eddy  Heat  Transport 

Tables  4.6A  and  4.6B  display  a  large  area  of  positive  correlations. 
Considering  all  the  latitudes  from  34°N  to  70®N,  we  see  that  the  area  is 
centered  at  zero  degrees  spatial  lag.  This  indicates  that  it  is  an  a  priori 
probability  peak,  and  the  contoured  area  is  significant  at  the  99%  confi¬ 
dence  level.  The  significance  of  the  correlations  is  strengthened,  since 
adjacent  latitudes  display  the  same  relationship.  The  analysis  discloses 
that  the  peak  transient  eddy  heat  transports  occur  at  approximately  the 
same  location  as  the  peak  baroclinic  instability.  This  observation  sup¬ 
ports  the  theory  that  transient  eddy  heat  transports  are  in  part  caused  by 
and  exist  to  alleviate  baroclinic  instability  (i.e.  mean  tropospheric  tem¬ 
perature  gradients) .  Obviously,  the  analysis  indicates  a  direct  relation¬ 
ship  between  these  two  quantities. 

At  latitudes  66°N  and  70°N,  a  negative  correlation  peak  exists  at 
approximately  80  degrees  spatial  lag.  This  observation  can  be  explained 
by  looking  over  the  whole  latitude  circle.  Clearly,  a  two-wave  pattern 
exists  for  both  quantities  at  these  latitudes,  and  this  observation  is 
supported  by  their  autocorrelations.  At  latitudes  34°N  and  38°N,  a  nega¬ 
tive  peak  occurs  between  290  and  310  degrees  spatial  lag.  Again,  this  can 
be  explained  by  the  cyclic  behavior  of  both  quantities.  In  this  case ,  both 
quantities  have  a  single  wave  pattern  around  the  latitude  circles.  Another 
interesting  observation  is  that  the  area  of  positive  correlation  appears 
to  have  a  slope  from  low  latitudes  with  a  small  positive  lag,  to  high  lati- 
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tudes  with  a  small  negative  lag.  This  apparent  slope  is  not  clearly  rep¬ 
resented  in  all  the  individual  monthly  correlation  arrays,  and  may  be  just 
a  result  of  the  averaging  process.  Obviously,  more  observations  would  be 
necessary  to  draw  any  conclusions  about  the  validity  of  this  apparent 
slope. 

The  individual  monthly  correlation  arrays  are  interesting  in  this 
case.  For  January  1973  and  1975,  the  peak  correlation  coefficients  occur¬ 
red  at  the  latitude  where  the  transient  eddy  heat  transport  had  its  largest 
value.  In  both  cases,  the  correlation  coefficient  was  greater  than  .70. 
This  observation  did  not  occur  in  January  1974.  In  fact,  this  month  was 
quite  different,  having  much  weaker  correlations  for  almost  all  latitudes. 
Perhaps  this  is  because  the  transient  eddy  heat  transports  were  much  weaker 
in  January  1974  (see  Table  2.4).  In  general,  the  correlation  coefficients 
appear  to  increase  when  the  transient  eddy  heat  transports  become  larger. 

4.7  Meridional  Surface  Temperature  Gradient  versus  Total  Meridional 


The  result  of  this  correlation,  displayed  in  Tables  4.7a  and  4.7B, 
is  one  we  might  have  expected.  The  correlation  array  appears  to  be  simi¬ 
lar  to  an  average  of  the  correlation  arrays  obtained  with  the  surface  tem¬ 
perature  gradient  versus  the  two  components  of  the  total  eddy  heat  trans¬ 
ports.  Since  the  stationary  component  displayed  no  clear  relationship  and 
the  transient  component  indicated  at  best  a  weak  relationship,  we  should 
have  expected  only  a  very  weak  relationship.  The  analysis  clearly  indi¬ 
cates  this  logic  to  be  true.  It  is  interesting  to  observe  that  no  differ¬ 
ent  relationships  are  suggested  when  we  combine  the  two  eddy  components. 
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We  do  have  a  few  peak  correlation  coefficients;  however,  they  reflect  the 
relationships  already  discussed.  The  relationship  between  the  surface 
temperature  gradient  and  the  transient  eddy  heat  transport  is  greatly  re¬ 
duced,  if  not  eliminated,  when  both  components  are  considered  together. 

This  indicates  that  the  only  relationship  that  exists  between  the  surface 
temperature  gradient  and  eddy  heat  transports  is  the  weak  relationship 
with  the  transient  component. 

4.8  Two-Layer  Model  Stability  Parameter  versus  Total  Meridional  Eddy 
Heat  Transport 

Tables  4.8a  and  4.8B  contain  the  correlation  analysis.  Again, 
there  is  no  different  relationship  suggested  by  this  analysis  other  than 
those  already  mentioned.  The  correlation  array  appears  to  be  like  an  aver¬ 
age  of  the  correlation  arrays  obtained  with  the  model  stability  parameter 
versus  the  two  components  of  the  total  eddy  heat  transport.  In  general, 
the  correlation  coefficients  are  smaller  than  those  obtained  with  the  model 
stability  parameter  versus  the  transient  eddy  heat  transport,  except  at 
low  latitudes  and  small  spatial  lags.  The  reason  that  the  correlation 
coefficients  are  not  reduced  in  this  location  is  because  of  the  similarity 
between  the  stationary  and  transient  components  (see  Table  4.2A).  Thus, 
the  analysis  indicates  that  the  only  clear  relationship  is  between  the 
model  stability  parameter  and  the  transient  eddy  heat  transport. 
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CHAPTER  FIVE 
SUMMARY  AND  CONCLUSIONS 

5.1  Assessment  of  Results 

In  Chapter  Four,  we  looked  for  possible  relationships  involving 
the  meridional  stationary  eddy  heat  transport.  The  analysis  indicated 
that  no  relationship  existed  between  the  meridional  stationary  eddy  heat 
transport  and  the  meridional  surface  temperature  gradient.  Also,  no  con¬ 
vincing  relationship  was  discovered  between  the  stationary  eddy  heat  trans¬ 
port  and  baroclinic  instability  (i.e.  the  two-layer  model  stability  para¬ 
meter)  .  rhis  indicates  that  the  position  and  strength  of  the  stationary 
eddy  heat  transports  are  primarily  caused  by  other  factors,  such  as  topo¬ 
graphy. 

The  meridional  transient  eddy  heat  transports,  on  the  other  hand, 
display  a  spatial  relationship  with  strong  meridional  surface  temperature 
gradients.  A  stronger  relationship  was  indicated  by  the  analysis  between 
the  transient  eddy  heat  transport  and  the  two-layer  model  stability  para¬ 
meter.  This  suggests  that  meridional  temperature  gradients  act  as  a  di¬ 
rect  forcing  agent  to  the  transient  eddy  heat  transports.  The  two-layer 
model  stability  parameter  has  a  stronger  relationship  with  the  transient 
eddy  heat  transport,  and  this  is  logical,  since  the  model  stability  para¬ 
meter  reflects  the  average  meridional  temperature  gradient  for  the  tropo¬ 
sphere.  The  relationship  between  the  model  stability  parameter  and  the 
transient  eddy  heat  transport  will  be  investigated  further  in  section  5.2. 

An  interesting,  but  not  totally  unexpected  observation  occurred 
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when  we  considered  the  total  eddy  heat  transport.  No  new  relationships 
were  suggested  by  the  analysis,  and  the  relationship  was  weaker  than  that 
between  the  transient  eddy  component  and  the  meridional  temperature  gra¬ 
dients.  The  relationship  was  weaker  because  of  the  dissimilarity  between 
the  stationary  and  transient  components.  The  dissimilarity  between  the 
two  eddy  components  suggests  that  they  are  primarily  caused  by  unrelated 
factors.  This  indicates  that  it  is  necessary  to  separate  the  total  eddy 
heat  transport  into  its  two  eddy  components,  when  investigating  for  possi 
ble  relationships. 

5.2  Determining  the  Two-Layer  Model  ..ability  Parameter  and  Meridional 
Transient  Eddy  Heat  Transport  Relationship 

•  The  transient  eddy  heat  transport  had  its  largest  values  at  lati¬ 
tudes  42°N  and  46°N  (see  Table  2.4).  looking  at  Table  4.6A,  the  largest 
correlation  coefficients  at  these  latitudes  occurred  with  a  spatial  lag 
of  ten  degrees.  For  these  reasons,  we  used  these  two  locations  to  inves¬ 
tigate  an  empirical  relationship  between  the  model  stability  parameter 
and  the  transient  eddy  heat  transport. 

5.2.1  Investigating  a  Linear  Relationship 

Figure  5.1  contains  a  graph  of  the  model  stability  parameter  ver¬ 
sus  the  transient  eddy  heat  transport  ac  latitude  42°N  and  with  a  spatial 
lag  of  ten  degrees.  In  this  graph,  all  three  months  of  data  were  used. 
This  provided  us  with  108  points  to  determine  a  linear  relationship.  A 
linear  regression  was  performed  in  order  to  determine  the  best  fit  line. 
The  equation  of  this  line  was 


i 
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Transient  Eddy  Heat  Transport  Cl  O'  *  calorie s/day/5*  longitude) 

Figure  5.1*  Plot  of  all  three  months  for  latitude  2*2*N  and  ten  degrees  spatial 
lag.  Each  circle  around  a  point  indicates  another  point  at  that 
location. 
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F  =  3.57 (MSP)  +  1.11 

17 

where  F  is  the  transient  eddy  heat  transport  in  units  of  10  calories 
per  day  per  five  degrees  longitude  and  (MSP)  is  the  model  stability  para¬ 
meter  in  units  of  meters  per  second.  looking  at  the  graph  in  Figure  5.1, 
we  see  a  large  scatter  around  the  line.  The  correlation  coefficient  be¬ 
tween  the  two  quantities  was  .42,  indicating  that  only  18%  of  the  varia¬ 
tion  of  the  transient  eddy  heat  transport  is  accounted  for  by  differences 
in  the  model  stability  parameter. 

The  same  calculation  was  performed  for  latitude  46°N  and  a  spatial 
lag  of  ten  degrees.  Figure  5.2  contains  the  graph  and  best  fit  line.  The 
equation  of  the  line  was 

F  *  4 (MSP)  +  1.48 

The  correlation  coefficient  in  this  case  was  .40;  therefore,  only  16%  of 
the  variation  of  the  trains ient  eddy  heat  transport  is  accounted  for  by  the 
differences  in  the  model  stability  parameter.  Obviously,  only  a  weak 
linear  relationship  is  indicated  in  both  cases.  The  accuracy  of  the  fit 
is  questionable  because  of  the  large  scatter  around  the  lines.  It  is  in¬ 
teresting  to  observe,  however,  that  the  equations  are  quite  similar  for 
the  two  different  latitudes. 

The  large  scatter  about  the  lines  may  be  due  at  least  in  part  to 
the  smoothing  of  the  model  stability  parameter.  If  we  had  not  smoothed 
the  model  stability  parameter  or  smoothed  this  parameter  in  a  different 
manner,  the  scatter  might  have  been  reduced.  On  the  other  hand ,  the  scatter 
may  well  be  due  to  other  important  factors  that  affect  the  transient  eddy 


Plot  of  all  three  months  for  la  itude  1|6*N  ancl  ten  degrees  spatial 
lag.  Each  circle  around  a  point  indicates  another  point  at  that 
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heat  transport.  Other  possible  factors  could  not  be  investigated  with  our 
present  data. 


5.2.2  Investigating  a  Power  Relationship 

The  scatter  diagrams  in  Figures  5.1  and  5.2  did  not  clearly  indi¬ 
cate  a  power  relationship;  however,  our  investigation  would  not  be  complete 
without  testing  this  possibility.  In  order  to  test  for  a  power  relation¬ 
ship,  we  added  Uc  to  the  model  stability  parameter.  This  made  all  these 
values  positive,  so  a  log  relationship  could  be  used.  The  critical  wind 
shear  is  essentially  a  constant  with  latitude  and  is  related  to  the  smooth¬ 
ing  distance,  Xq,  by  Stone  (1978) : 


U 

c 


,  flaTT2 
10368 


cos3  |J> 


(5.1) 


where  ft  is  the  earth's  angular  velocity,  a  is  the  mean  radius  of  the  earth, 

$  is  the  latitude,  and  X  is  the  smoothing  distance  in  increments  of  five 

o 

degrees  longitude.  Using  this  equation,  Uc  was  calculated  as  11.6  meters 
per  second  and  9.5  meters  per  second  for  latitudes  42°N  and  46°N,  respec¬ 
tively.  The  relationship  we  will  be  testing  is 

F  =  b (MSP+U  )d  (5.2) 

c 

where  F  is  the  transient  eddy  heat  transport,  MSP  is  the  model  stability 
parameter,  Uc  is  the  critical  wind  shear,  and  b  and  d  are  constants  to  be 
determined.  Taking  the  log  of  equation  5.2,  we  have 

log  F  *  d  log  (MSP+Uc)  +  log  b  (5.3) 


*\ 

k 
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Equation  5.3  is  in  the  form  of  a  line,  and  we  can  perform  a  linear  regres¬ 
sion  on  log  F  and  log  (MSP+U  )  to  calculate  the  values  of  b  and  d.  All  of 

c 

the  data  points  were  not  used  in  this  investigation,  because  zero  or  nega¬ 
tive  values  of  the  transient  eddy  heat  transport  could  not  be  employed  in 
equation  5.3. 

Figure  5.3  displays  a  graph  of  (MSP+lM  versus  F  for  our  three 
months  plotted  on  log  x  log  paper  for  latitude  42 °N  and  a  spatial  lag  of 
ten  degrees.  Performing  a  linear  regression  on  the  data,  two  lines  were 
calculated.  Line  1  is  the  best  fit  line  calculated  by  minimizing  the 
variation  of  the  transient  eddy  heat  transport  about  a  line,  and  has  the 
equatior 

86 

F  =  . 59 (MSP+U  )‘ 
c 

As  before,  (MSP+Uc)  is  in  units  of  meters  per  second  and  F  is  in  units  of 
17 

10  calories  per  day  per  five  degrees  longitude.  Line  2  is  the  best  fit 
line  determined  by  minimizing  the  variation  of  (MSP+Uc)  about  a  line,  and 
its  equation  is 

F  =  (2.6xl0-9) (MSP+U  )8'3 

c 

If  there  were  an  exact  relationship  between  these  two  quantities,  it  would 
lie  between  or  on  one  of  these  lines.  The  large  difference  between  these 
two  lines  is  due  to  the  large  scatter.  The  correlation  coefficient  be¬ 
tween  log  F  and  log  (MSP+Uc)  is  .32,  indicating  that  only  10%  of 

the  variation  of  log  F  is  explained  by  the  differences  in  log  (MSP+U  ) . 

c 


Transient  Eddy  Heat  Transport  (101 '  calories/day/5“  longitude) 
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A  test  of  the  scatter  about  each  line  indicated  that  line  2  was  by  far  the 
best  fit  to  our  data. 

Figure  5.4  is  a  graph  of  (MSP+U^)  versus  F  at  latitude  46°N  and  a 
spatial  lag  of  ten  degrees.  Line  1  was  calculated  by  minimizing  F  about  a 
line,  and  has  the  equation 

F  =  .86  (MSP+U  ) *91 
c 

Line  2's  equation  is 

F  =  (4.29xl0-5) (MSP+U  )5,1 

c 

and  was  calculated  by  minimizing  (MSP+U^)  about  a  line.  Again,  a  scatter 

test  was  performed,  and  strongly  indicated  that  line  2  was  the  best  fit  to 

our  data-  The  correlation  coefficient  for  log  F  and  log  (MSP+U  )  was  .42; 

c 

therefore,  18%  of  the  variation  of  log  F  is  explained  by  differences  in 
log  (MSP+U^) .  Since  the  correlation  between  these  quantities  is  larger 
than  for  latitude  42°n,  the  uncertainty  in  the  relationship  and  the  differ¬ 
ence  between  the  two  equations  is  smaller.  The  equation  for  line  2  above  had 
the  best  fit  of  the  relationships  determined  and  agrees  with  Held's  (1978)  sug¬ 
gested  fifth  power  dependence.  Still,  as  in  the  linear  investigation,  the 
data  is  not  conclusive  enough  to  give  a  specific  relationship. 

5.3  Areas  for  Further  Investigation 

As  we  indicated  in  section  5.2.1,  our  smoothing  of  the  model  sta¬ 
bility  parameter  may  have  caused  some  of  the  scatter  in  our  graphs.  In 
smoothing,  we  may  have  eliminated  some  of  the  variation  of  the  model  sta¬ 
bility  parameter  related  to  the  variation  of  the  transient  eddy  heat  trans¬ 
port,  thus  decreasing  their  correlation.  An  investigation  with  an  un¬ 
smoothed  model  stability  parameter  or  the  model  stability  parameter  smoothed 
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Transient  Eddy  Heat  Transport  (10  '  calories/day/5*  longitude) 
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MSP  +  Uc  (me ters/sec) 


Figure  5 *4*  Log  base  ten  plot  of  all  three  months  for 
latitude  14.6*  li  and  ten  degrees  spatial  lag. 
Each  circle  around  a  point  indicates 
another  point  at  that  location. 
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in  a  different  manner  may  prove  more  fruitful.  Another  possibility  is 
that  the  relationship  may  be  improved  by  performing  the  same  smoothing 
on  both  quantities. 

Also,  the  data  should  be  investigated  over  a  much  longer  time 
period.  Our  observations  could  then  be  checked  to  see  if  they  persist, 
and  interannual  variations  could  be  investigated.  With  a  larger  data  base, 
the  a  posteriori  correlation  peak  between  the  surface  temperature  gradient 
and  the  transient  eddy  heat  transport  could  be  investigated  more  completely. 
Different  months  should  also  be  studied,  to  see  if  the  relationships 
change  with  season. 

” ^rthermore ,  a  time  series  of  the  local  model  stability  parameter 
and  the  transient  local  eddy  heat  transport  should  be  examined,  so  the 
importance  of  time  lags  between  the  two  quantities  could  be  investigated. 

We  attempted  to  do  this  over  a  ten  day  period  in  January  1973,  with  oui 
quantities  calculated  every  twelve  hours.  Since  the  raw  data  was  only 
measured  every  twelve  hours,  we  could  not  separate  the  eddy  components. 

This  only  allowed  us  to  compare  the  model  stability  parameter  with  the 
total  eddy  heat  transport.  The  only  peak  in  the  correlation  analysis  oc¬ 
curred  in  mid-latitudes  at  small  spatial  lags,  similar  to  what  we  discov¬ 
ered  in  our  previous  investigation  of  these  two  quantities.  The  peak  was 
always  pr°sent  in  this  location,  regardless  of  the  time  lag,  so  the  analy¬ 
sis  provided  no  new  insights.  This  result  may  have  been  caused  by  the 
dampening  effect  of  the  stationary  eddy  component  on  the  relationship  be¬ 
tween  the  transient  eddy  heat  transport  and  the  model  stability  parameter, 
when  both  eddy  components  are  added  together.  If  we  had  considered  a 
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longer  time  period  and  calculated  the  transient  eddy  heat  transport  for 
every  two  or  three  days,  then  performed  our  analysis  between  the  model 
stability  parameter  and  transient  eddy  heat  transport,  we  may  have  dis¬ 
covered  a  greater  insight  in  their  relationship. 

Finally,  relationships  between  meridional  eddy  heat  transports 
and  other  factors  should  be  considered.  In  our  investigation,  the  corre¬ 
lation  between  the  model  stability  parameter  and  the  transient  eddy  heat 
transport  was  weak,  although  significant.  Possibly  meridional  temperature 
gradients  are  not  the  only  important  forcing  agents  for  meridional  trans¬ 
ient  eddy  heat  transports.  Other  factors  that  may  play  an  imortant  role 
in  the  f -^mation  and  maintenance  of  transient  eddy  heat  transports  are 
topography,  latitudinal  variations  of  temperature,  etc. 
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